Dynamics of Microtubule Instabilities 
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We investigate an idealized model of microtubule dynamics that involves: (i) attachment of 
guanosine triphosphate (GTP) at rate A, (ii) conversion of GTP to guanosine diphosphate (GDP) 
at rate 1, and (iii) detachment of GDP at rate /x. As a function of these rates, a microtubule can 
grow steadily or its length can fluctuate wildly. For fj, — 0, we find the exact tubule and GTP cap 
length distributions, and power-law length distributions of GTP and GDP islands. For /i = oo, we 
argue that the time between catastrophes, where the microtubule shrinks to zero length, scales as 
e A . We also discuss the nature of the phase boundary between a growing and shrinking microtubule. 

PACS numbers: 87.16.Ka, 87.17.Aa, 02.50.Ey, 05.40.-a 



Microtubules are linear polymers of the protein tubu- 
lin that perform major organizational tasks in living cells 
[l|, Q • They provide transport tracks for molecular ma- 
chines [H, H|, ano - m ove cellular structures during cellular 
processes such as reproduction [2j,|5|. A surprising feature 
of microtubules is that they remain out of equilibrium un- 
der fixed external conditions and can undergo alternating 
periods of growth and rapid shrinking (Fig. [1]) . 
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FIG. 1: Schematic illustration of the growth and catastrophic 
shrinkage of a microtubule as a function of time (adapted 
from 

These sudden polymerization changes are driven by the 
interplay between growth by the attachment of guanosine 
triphosphate tubulin complexes (GTP) at the free end 
[J,IH @, 0] i irreversible hydrolysis of GTP into guanosine 
diphosphate (GDP) anywhere in the tubule, and subse- 
quent GDP detachment from the free end of the tubule. 
One avenue of theoretical work on this dynamical insta- 
bility is based on detailed models of mechanical stabil- 
ity [1, Q. For example, a detailed stochastic model of 
a microtubule that includes all the thirteen constituent 
protofilaments has been investigated in Ref. 0]. By us- 
ing model parameters that were inferred from equilibrium 
statistical physics, VanBuren et al. [§] found character- 
istics of the tip evolution of a microtubule that agreed 
with experimental data [lfj. 

Another approach for modeling the dynamics of mi- 
crotubules is based on effective two-state models that 
describe the dynamics in terms of a switching between 
a growing and a shrinking state [Til. [l2l [Tl. fl4 [Tol [T(| . 
The essence of many of these models is that a microtubule 



exists either in a growing phase (where a GTP cap exists 
at the end of the microtubule) or a shrinking phase (with- 
out a GTP cap), and that there are stochastic transitions 
between these two states. By tuning parameters appro- 
priately, it is possible to reproduce the phase changes be- 
tween the growing and shrinking phases of microtubules 
that have been observed experimentally J3J . In another 
important contribution, Flyvbjerg et al. [171 ] constructed 
an effective continuous theory to describe the dynamics 
of the cap length. 

What is missing in these models is a clear connection 
between microscopic parameters and the evolution of the 
spatial structure of the microtubule. Most of the models 
discussed thus far are too complicated to permit a com- 
plete solution. Thus the present work is motivated by 
the goal of formulating and solving a minimal and ideal- 
ized coarse-grained model of microtubule dynamics that 
incorporates the main features of growth and shrinking. 
We view this idealized model as an Ising-like descrip- 
tion that captures the most interesting features of mi- 
crotubule dynamics even though the connection between 
model parameters and experimental variables is indirect. 
One of the main advantages of the model is its simplicity 
so that many of its geometrical and dynamical features 
can determined analytically. These solutions reveal many 
rich geometrical and time-dependent features that should 
help our understanding of real microtubule dynamics. 

The model studied in this paper [l8[ treats a micro- 
tubule as a linear polymer that consists of GTP or GDP 
monomers that we denote as + and — , respectively. To 
emphasize this connection between chemistry and the 
model, we write the former as GTP+ and the latter as 
GDP - . The state of a microtubule evolves by the follow- 
ing steps: 

1. Attachment of GTP+ at the end of a microtubule: 

!•••) =H h) rate A 

2. Independent conversion of each GTP + to GDP - : 

!••• + •••) =► | ) rate 1. 

3. Detachment of a GDP - from the microtubule end: 

| • • • — ) =>■ | ■ ■ • ) rate /i. 
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Here the symbols | and ) denote the terminal and the 
free end of the microtubule. It is worth mentioning that 
these steps are similar to those in a recently-introduced 
model of DNA sequence evolution [l9[ , and that some of 
the results about the structure of DNA sequences seem 
to be related to our results about island size distributions 
in microtubules. 

With the above steps for the evolution of a micro- 
tubule, we find that the (A, fi) phase plane separates into 
a region where the microtubule grows, on average, and a 
phase where the microtubule length is finite. On the 
boundary between these two regions, the microtubule 
length fluctuates wildly. To understand these different 
phases, we first focus on the extreme cases of no detach- 
ment fi — and infinite detachment rate fi = oo, where 
we can analytically solve the microtubule structure. 

No Detachment : Here GTP + monomers attach at rate 
A and convert to GDP - at rate 1. The probability H^{t) 
that the microtubule consists of A GTP+ monomers at 
time t evolves according to 
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The loss term (A + \)Hn accounts for the attachment of 
a GTP+ at the end of the microtubule of length A with 
rate A and the conversion events GTP + — > GDP - that 
occur with total rate A. The gain terms can be explained 
similarly. This equation can be solved by the generating 
function method and the final result is 
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From this Poisson distribution, the mean number of 
GTP + monomers and its variance are 



(N) = (A 2 ) - (AO 2 = A(l -e -t ). 



(3) 



Similarly, the length distribution P(L,i) of the micro- 
tubule evolves according to the master equation 



dP{L,t) 
dt 



\[P{L - l,t) - P(L,t)] 



(4) 



For the initial condition P(L,0) — <5_l.o, the solution is 
again the Poisson distribution 



P(L,t) 
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(5) 



from which the growth rate of the microtubule and the 
diffusion coefficient of the tip are V = A and D = A/2. 

Because of the conversion of GTP + to GDP - , the tip 
of the microtubule is comprised predominantly of GTP + , 
while the end exclusively consists of GDP - . The region 
from the tip until the first GDP - is known as the cap 
(Fig. [5]) and it plays a fundamental role in microtubule 
function. We use a master equation approach to deter- 
mine the cap length distribution (20| . 



Consider a cap of length k. Its length increases by 1 
due to the attachment of a GTP+ at rate A. The con- 
version of any GTP+ into a GDP - at rate 1 reduces the 
cap length from k to an arbitrary value s < k. These 
processes lead to the following master equation for the 
probability n k that the cap length equals k: 



hk = A(rifc_ 



rik) 



krik + ^ 



(6) 
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Equation © remains valid for k = if we set n_i = 0. 
We solve for the stationary distribution by summing the 
first k — 1 of Eqs. © with set to zero to obtain 



rik-i = 



E 1 

s>k 



(7) 



The cumulative distribution, Nk = X) s >fc ris ' thus satis- 
fies the recursion Nk — XNk-i/(k + A). Using the nor- 
malization Ao = 1 and iterating, we obtain the solution 
in terms of the Gamma function [21(: 
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r(fc + 1 + A) 

Hence the cap length distribution is 

r(i+A) 



T(fc+2 + A) 



(fc + l)A fe 



(8) 



(9) 



From this distribution, we find that in the realistic 
limit of A ^> I the average cap length is given by 



(k) 



as A 



oo 



(10) 



Thus even though the average number of GTP + 
monomers equals A, only -\/A of them organize themselves 
into the microtubule cap. 
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FIG. 2: Representative configuration of a microtubule, with 
a cap of length 4, three GTP + islands of lengths 1, 3, and 2, 
and three GDP - islands of lengths 3, 2, and a "tail" of length 
5. The rest of the microtubule consists of GDP - . 

At a finer level of resolution, we determine the distri- 
bution of island sizes in the GTP + -populated zone of a 
microtubule (Fig. [J). When A — > oo, both the length of 
the cap and the length of the region that contains GTP+ 
become large and a continuum description is appropri- 
ate. Since the residence time of each monomer increases 
linearly with distance from the tip and the conversion 
GTP + — * GDP - occurs independently, the probability 
that a monomer remains a GTP + decays exponentially 
with distance from the tip. This fact alone is sufficient 
to derive all the island distributions. 
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Consider first the length £ of the populated zone 
(Fig. [5]). For a monomer at a distance x from the tip, its 
residence time is r = x/X for large A. Thus the probabil- 
ity that this monomer remains a GTP+ is e~ T = e~ x / x . 
We thus estimate I from the extremal criterion [221 



l = ^e- K / A = (l-e- 1 / A )- 1 £ 



-e/x 



(11) 



x>£ 



that merely states that there is of the order of a single 
GTP + further than a distance £ from the tip. When A 
is large, (1 — e -1 / A ) -1 — > A, and the length of the active 
region scales as 

£ = AlnA (12) 
The probability that the cap has length k is given by 



(1 _ e -(fc+i)A) JJ e -o/x, 

3 = 1 

The product ensures that all monomers between the tip 
and a distance k from the tip are GTP + , while the pref- 
actor gives the probability that a monomer a distance 
k + 1 from the tip is a GDP - . Expanding the prefactor 
for large A and rewriting the product as the sum in the 
exponent, we obtain 



1 



-fe(fe+l)/2A 



A 



(13) 



a result that agrees with the large-A limit of the exact 
result for n k in Eq. Q. 

Similarly, the probability to find a positive GTP + is- 
land of length k that occupies sites x + l,x + 2, . . . ,x + k 
is 

k 

(1 _ e -x A)(i - e -(*+fc+i)A) JJ e -0+j)A_ (14) 

j'=i 

The two prefactors ensure that sites x and x + k + 1 are 
GDP - , while the product ensures that all sites between 
x + 1 and x + k are GTP+. 

Most islands are far from the tip and they are short, 
k < x, so that (HU) simplifies to (I- e- x / x ) 2 e~ kx / x . The 
total number of GTP+ islands of length k is obtained by 
summing this island density over all x. Since A > 1, we 
replace the summation by integration and obtain 



J* 



dx (1 - e ^A)2 e -fcz/A 



2A 



fc(fc+l)(fc+2) 



(15) 



By similar reasoning, we find that the density of neg- 
ative GDP - islands of length k <C x with one end at 
x is given by e -2x / A (l — e -x A) fc . xhe total number of 
negative islands of length k is therefore 

poo \ 

* = /. ^e-^(l-e-A)- = (t + 1)(t + 2) .(16) 



Again, we find a power-law tail for the GDP - island size 
distribution, but one that is much broader than the cor- 
responding GTP+ distribution. Strikingly, these two dis- 
tributions are of the same form as those found for the de- 
gree distribution of growing networks (23l . l2~j ] . The total 
number of GDP - monomers within the populated zone 
is Ylk>i kJk- While this sum formally diverges, we in- 
voke the upper size cutoff, /c* ~ A that again follows from 
an extremal criterion [22j | to obtain ^2 k >ikJ k — A In A. 
Since the length of the populated zone £ ~ A In A, we 
thus see that this zone predominantly consists of GDP - 
islands. 
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FIG. 3: Cartoon of a microtubule. The GTP + regions (dark) 
get shorter further from the tip that advances as Xt, while the 
GDP - regions (light) get longer. 

In analogy with the cap, consider now the "tail" — the 
last island of GDP - within the populated zone (Fig. [2]). 
The probability m k that it has length k is 



m k 



= e-"/ x {l-e- l l x ) k . 



(17) 



Hence, by summing the geometric series and using I = 
A In A from (I12p , the average length of the tail is 

(k) = km k = e e/x - 1 — > A. (18) 
fc>i 

Thus the tail is (on average) much longer than the cap. 
In summary, the microtubule consists of GTP + islands 
that systematically get shorter away from the tip and 
vice versa for GDP - islands (Fig. [3]). 

Instantaneous Detachment : For fj, > 0, a micro- 
tubule can recede if the monomer(s) at its tip are GDP - . 
If there are many such GDP - end monomers the micro- 
tubule length can fluctuate wildly under steady condi- 
tions. Here we focus on the limiting case of instantaneous 
detachment, fi = oo. As soon as the monomer at the tip 
changes from a GTP + to a GDP - , this monomer and 
any contiguous GDP - monomers detach immediately. 
We term such an event an avalanche of size k. If the 
avalanche encompasses the entire microtubule, we have a 
global catastrophe. The probability for such a catastrophe 
to occur is 



e(A) 



l + A 



71=1 



-n/X\ 



(19) 



The factor (1 + A) -1 gives the probability that the 
monomer at the tip converts to a GDP - before the next 
attachment event, while the product gives the probability 
that all other monomers in the microtubule are GDP - . 
In principle, the upper limit in the product is set by the 
microtubule length. However, for n > A, each factor in 
the product is close to 1 and the error made in extending 
the product to infinity is small. 
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The asymptotic behavior of the infinite product in 
Eq. (|19[) is found by first expressing it in terms of the 
Dedekind rj function [25| 



T](z) = e™ z/12 Y[{1 - e 27rmz ), 



(20) 



n=l 



and then recalling a remarkable identity satisfied by this 
function, rj(— 1/z) = \ / —izr](z). Using these facts, the 
probability of a catastrophe is given by 



e(A) 



l + A 



r 2 A/6 e l/24A "Q ^ _ e -47T 2 An) 



n>l 



27[ e ^ 2 A/6 



(21) 



Since the time between catastrophes scales as the inverse 
of the occurrence probability, this inter-event time be- 
comes astronomically long for large A. 




FIG. 4: Schematic phase diagram for a microtubule in the 
A-/j parameter space. 

Similar arguments can be used to calculate the size 
distribution of finite avalanches. We first notice that the 
configuration | 4 — • ■ ■ — h) occurs with probability 

k-l 

k-l 



no 



-n/X\ 



71=1 



Then multiplying by the probability that the monomer at 
the tip converts before the next attachment event gives 
the probability for an avalanche of size k: 



k-l 



Ai 



(l + A)- 1 JJ(l- e -»A) (22) 
Using 1 - e~ n / x « n/A - n 2 /(2A 2 ), Eq. ([22]) becomes 



k-l 



A k = \- k r(k) n (i 



2A 



\- k T{k)e 



-k 2 /4\ 



General Growth Conditions: For arbitrary attach- 
ment and detachment rates A and /i, we can approxi- 
mately map out the boundary between different phases of 



microtubule dynamics in the A, \x parameter space. For 
A, [i <C 1, the unit conversion rate GTP+ — > GDP~ is 
much faster than the rates A, fi of the other dynamical 
microtubule processes. Hence we assume that conversion 
is instantaneous. Consequently, a microtubule consists 

of a string of GDP~ monomers | • • ■ ) in which the 

tip advances at rate A and retreats at rate /i. Thus the 
tip performs a biased random walk, and its velocity is 
V(X, fi) = A— /i when A > fJ,. The boundary between the 
growing phase, where the average tubule length grows 
as Vt, and the compact phase, where the average tubule 
length is finite, is found by setting V = 0. This condition 
gives /i* = A when A <C 1. On this phase boundary, the 
average tubule length grows as \fi. 

For the phase boundary for large /i, consider first 
fi = oo. Since the leading behavior for the probability 
of a catastrophe scales as e _Tr A / 6 , the typical time be- 
tween catastrophes is e w A / 6 . Since V « A for large A, the 
typical length of a microtubule just before a catastrophe 
is again (ignoring all power-law factors) e* A / 6 . Suppose 
now the detachment rate /i is very large but finite. The 
microtubule is compact if the time to shrink a micro- 
tubule of length e" A / 6 //i is smaller than the time A -1 

required to generate a GTP+ by | ■ ) ==> | ■ • • — h) 

and thereby stop the shrinking. We estimate the loca- 
tion of the phase boundary by equating these two times 
to give fi* - e * 2x /s when A > 1 (Fig.|4|). 

To summarize, our model predicts rich microtubule dy- 
namics as a function of GTP + attachment and GDP~ 
detachment. In the growing phase, GTP + and GPD~ 
organize into alternating domains, with gradually longer 
GTP + domains and gradually shorter GDP~ domains 
toward the tip of the microtubule. The size distributions 
of these two species are exact power laws with respective 
exponents of 3 and 2. Between the limiting cases of a 
finite-length and a growing microtubule, its length can 
fluctuate wildly under steady external conditions. This 
unusual behavior emerges naturally in our model. From 
a simple probabilistic approach and in the limit instan- 
taneous detachment of GDP~ (/j = oo), we found that 
the time between catastrophes, where the microtubule 
shrinks to zero length, scales exponentially with the at- 
tachment rate A. Thus for large A, the microtubule will 
grow essentially freely for a very long time before under- 
going a catastrophe. 
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